Bonn TK 97-08 



QCD at large and short distances^ 



Vladimir Gribov 

Landau Institute for Theoretical Physics, Moscow 

and 

Research Institute for Particle and Nuclear Physics, Budapest 

and 

Institut fur Theoretische Kernphysik der Universitdt Bonn 

July 28, 1997 



Abstract 

The formulation of QCD which contains no divergences and no renor- 
malization procedure is presented. It contains both perturbative and non- 
perturbative phenomena. It is shown that, due to its asymptotically free 
nature, the theory is not defined uniquely. The chiral symmetry breaking 
and the nature of the octet of pseudo-scalar particles as quasi-Goldstone 
states are analysed in the theory with massless and massive quarks. The 
U(l) problem is discussed. 



1 Introduction 



In this paper I show how to formulate an asymptotically free theory in such a way 
that it includes perturbative and non-perturbative phenomena simultaneously. 

The idea is the following. Contrary to an infrared free theory, in an asymp- 
totically free theory the divergences prevent us from writing even perturbation 
expansions in a unique, well-defined way. We can, however, make use of the fact 
that divergences in the theory enter only the Green's functions and the vertices. 
On the other hand, knowing the Green's functions and the vertices we can express 
all the other amplitudes through them perturbatively in a unique way. Thus, we 
have to formulate equations for Green's functions and vertices in a form which 
does not contain any divergences. If this is done, the solutions of these equations 
will contain both perturbative and non-perturbative phenomena. 

*This work was completed in Germany under the Humboldt Research Award Program 
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To avoid technical complications, in the first section of the paper we will derive 
these equations in an Abelian theory in which usual perturbation theory is also 
well defined. In the second section we generalize the equations for a non- Abelian 
theory and discuss the connection between perturbative and non-perturbative 
effects. The equations have an integro-differential structure in which the asymp- 
totic behaviour of the Green's functions is defined by the boundary conditions. 
The main conclusion in this section is that in the region of large momenta the 
Green's functions of quarks and gluons contain additional parameters compared 
to normal perturbation theory which can be associated with different types of 
" condensates" . These non-perturbative parameters can be defined by solving the 
equations and finding non-singular solutions in the infrared region. A priori it is 
not clear what type of additional conditions have to be imposed on this system of 
equations in order to fix a non-singular solution. It can be, for example, the con- 
servation of the axial current or of other currents which is formally satisfied from 
the point of view of the Lagrangian but is not ensured because of the divergences. 

In the third section we will show that these equations make it possible, for 
the first time, to analyse the problem of spontaneous symmetry breaking in an 
asymptotically free theory and to find approximately the value of the critical 
coupling at which symmetry breaking occures. Also, they allow us to answer the 
fundamental question for asymptotically free theories, namely: how the bound 
states - the hadrons - have to be treated in these theories and how these bound 
states influence the equations for the Green's functions. 

The analysis leads to the following conclusion. The conditions for axial cur- 
rent conservation of flavour non-singlet currents (in the limit of zero bare quark 
masses) require that eight Goldstone bosons (the pseudo-scalar octet) have to be 
regarded as elementary objects with couplings defined by Ward identities. This 
is so in spite of the fact that the couplings of these states to fermions decrease at 
large fermion virtualities. 

The same analysis provides a new possibility for the solution of the U(l) 
problem. In this solution the flavour singlet pseudoscalar rj' is a normal bound 
state of qq without a point-like structure. The mass of this bound state is different 
from zero and can be calculated in the limit of massless quarks. For massive 
quarks the pseudoscalar octet becomes massive. Their masses, however, are not 
calculable in terms of bare quark masses because of logarithmic divergences and 
have to be regarded as unknown parameters which in the real case of confined 
quarks are defined by the self-consistence of the solution of the infrared problem. 
These states have an essential influence on the equations for the Green's functions, 
which, as it will be shown in the next paper, can be used constructively in solving 
the confinement problem if the effective coupling in the infrared region is not too 
large. In this case the integro-differential equations can be reduced to a system 
of non-linear differential equations and the theory looks like a theory of particle 
propagation in self-consistent fields defined by the Green's functions themselves 
(as it is the case in Landau's Fermi liquid theory). The self-consistent fields are 
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fields of gluons and 7r-mesons. 



2 Equations for Green's functions in QED 

In QED we have two Green's functions: that of the photon D^ u (k) and of the elec- 
tron G(q), and a vertex function T^(k, q). The photon Green's function is defined 
by the vacuum polarization operator Xl^ik) which can be expressed symbolically 
by a sum of Feynman diagrams 

IW*;) = eg{ 7 „C_> lu + 1»<S> 7, + • • '} (1) 

In order to obtain series not containing divergences we will try to consider the 
derivatives of H til/ (k) as functions of momenta. Due to current conservation, we 
have 

IW*) = {S^k 2 - k^K)U(k 2 ) (2) 

where U(k 2 ) contains only logarithmic divergences. 
Differentiating (Q) twice we will have 

d 2 U, v (k) = 6U(k 2 ) + (s^ - + 6)d^U(k 2 ); (3) 

we here use 

The second term in (|3|) does not contain divergences, the first one does. In order 
to obtain a finite expression, consider 

d^U^k) = -36^U(k 2 ) - 3^d^U(k 2 ). (4) 
Because of this, the quantity 

d 2 U, u (k)+2d,d a U au (k) = -6^d ? n(A; 2 )+ (s, v - (ds + 6)dtll(k 2 ) (5) 

must not contain divergences. 

In Feynman gauge D^ v {k 2 ) is equal 

n 1 _ 1 e 2 (k 2 ) 

M J ~ife2(l-n(Jfc2)) " k*el{k 2 ) { ' 

The diagrams contain only the product e^D^k) and therefore 

e 2 D, u (k) = ^e 2 (k 2 ) 
d,U(k 2 ) = -d,f 2 . (7) 



In first order this means 




the integrals on the right-hand-side are convergent. This fact has to hold in any 
order. Hence, we can include in (|8]) the exact electron Green's functions and the 
exact vertex functions and add all the corresponding diagrams which were not 
included. As a result, we can write 




where J2 P denotes the sum over the permutations of the indices similarly to @; 
a quantity p corresponds to each photon line. (^) is an equation for e 2 (in the 
form of series in e 2 ) provided that G(k) and T^(k, q) are known. 

In order to obtain equations for electron Green's functions and vertices we 
have to remember that these function can change rapidly even if e 2 is small be- 
cause they can have infrared singularities of the type e 2 In 9 ~" or ultraviolet 

singularities of the type (^) 7 - It is proven to be possible to arrange the differ- 
entiation in such a way that there will be an expansion only in e 2 . 

To write an equation for the fermion Green's function not containing any 
divergences, we have to differentiate twice the self-energy of the fermion or its 
inverse Green's function. Let us consider d flti G~ l (q)] the diagrammatic expression 
for G^ 1 will be the following. The simplest diagram is 



Diagrams of the next order are 



k 

>i — - — — — : v r- _|_ 



q-k'-k 



It can be easily shown that 



9 2 ^1- = -An 2 iS 4 (k) (10) 
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Using this equality, we have 

e 2 



d \/\/ = -g • — <■ — • — - -g < — • — ■ — • + 5 2 , 



k 2 



"7m- 



Restricting ourselves to <9 2 p-, we can write 

d 2 ai = -gol^Gol^, 

d 2 a 2 = -goi^God^ + god^Goin - go^G^ - ga^Goi^ 
and, consequently, 

d 2 G l = gd^Gd^G- 1 , Gf 1 + G 2 1 = G" 1 . (11) 

The term 5 2 has a contribution to 7^ \ " \ not containing overlapping diver- 
gences and is of the form 

_2fc CT _2K 

la la 

_ 2fc CT _ 2k' a 

la la 

For the calculation of higher order diagrams it is convenient to adopt the following 
principle. Beginning from the first point of interaction we shall relate the external 
momentum to the photon line: 



(13) 




q—q- 
q 



The next photon interaction can be expressed as 



now we relate the external momentum to the positron line. If the positron is 
emitting a photon, we keep sending the external momentum along the positron 
line up to its annihilation. As a result, we get a tree structure 




in which the initial photon line can end only at the final electron. It can be 
shown that every line in the diagram will be passed only once if the photons 
included in the fermion self energy are not taken into account. This means 
that in this approach exact electron Green's functions have to be used with bare 
vertices since the idea is basically the exclusion of overlapping divergences. Taking 
the second derivative of one of the photon propagators and restricting ourselves 
to the contribution 47r 4 5 4 (fc)g(0) we obtain both on the left-hand side and the 
right-hand side photon emission amplitudes with zero momentum. Due to gauge 
invariance, however, a zero momentum photon cannot change the state of the 
system (it is emitted from the external end). The photon emission amplitude 
equals 0) = d^G~ l (q). The final contribution of the differentiation has the 
structure flTTl) . The above statement, which essentially means that the emission 
of a zero momentum photon is determined by the total charge, can be proven by 
the Ward identity. 
Thus, we have 

d 2 G~ l = g^d.G^Gd.G- 1 + ( 16 ) 

The remaining diagrams contain first derivatives of photon and positron lines 
and second derivatives of positron lines. All these diagrams can be expressed 
in terms of the exact T^, G and D functions. They do not contain divergences 
except those graphs which correspond to the photon self energy. 
The first term in fll6| ) has the structure 

~g(0)M w (q,k = 0) 

where M uu (q, k = 0) is a quantity close to the Compton scattering amplitude in 
the sense that they would be equal if we differentiated all photon propagators 
including those inside the electron Green's function. (The factor | is due to the 
fact that the Compton amplitude contains the sum of diagrams with momenta 
k and — k). The Compton scattering amplitude M ul/ (q, k = 0) satisfies the Ward 
identity 

M w {q) = G- 1 d 2 GG- 1 = 2d IM G~ 1 Gd IM G- 1 = d 2 G~ x 
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which has the simple diagrammatical meaning 




The second term on the right-hand side of ([LTD represents the contribution of 
photons which, if we carried out the mentioned differentiation, would correspond 
to photons inside the Green's function and which are not present in M vv . Hence, 

M w (q, 0) = 2d tl G~ 1 Gd ll G- 1 . 

The remaining diagrams in ( |T6| ) contain first derivatives of photon and posi- 
tron lines and second derivatives of positron lines. All these diagrams can be 
expressed in terms of the exact T M , G and D functions. They do not contain 
divergences except those graphs which correspond to the photon self energy. All 
these diagrams are of the form 




and correspond to the photon line with the exact Green's function (J^) which 
equals 47r 2 g(/c)p-. By differentiating the photon lines we have calculated the 
contribution of d 2 -^. So, there remains 

1 = I Q^Wk) - 9(0))^ (19) 

where we have introduced g(0) in order to avoid a contribution from 5 4 (k). This 
expression contains logarithmic corrections coming from the ultraviolet region 
{k > q). Replacing g(k) — g(0) by g(q) — g(0) + g(k) — g(q) and performing an 
integration by parts, we get 

l=-(g(q)-g(0))M„v(q,0)+ 
+ J^(9(k)-g(q))^d 2 M vv (q,k). (20) 

The first term in ( f20|) can be explicitly calculated while the second one does not 
contain any logarithms because of the presence of the difference g(k) — g(q). We 
can rewrite (|20D in the form 

X = (g(q) - gmd.G-'Gd.G- 1 + 6 1 ; 
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consequently, d 2 G 1 can be expressed as 



d 2 G-\q) = g(q)d fl G- 1 Gd,G- 1 + 6 1 + 5 2 



(21) 



where 5i is defined by (^) and <5 2 contains first order derivatives of the photon 
and positron lines and second order derivatives of the positron lines, as it has been 
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explained above. The first term contains all singularities of the types a In q ~" 
and (^Y. 

3 Equations for the vertex function and for the 
amplitudes of interaction with the external 
field 

To obtain the equation for the vertex function, let us express T^(p, q) as a set of 
diagrams containing exact Green's functions: 





+ 



(22) 



We shall relate the external momenta to the lines in the diagrams in the same 
way as we did when we derived the equation for the Green's function. Calculating 
the second order derivative in q we obtain 

d 2 q T,(p, q) = g(0)MZ(p, q, k)\k=o + ■ - ■ (23) 

where M„ u (p, q, k)\k=o is the contribution to the amplitude of the process 




which corresponds to the emission of photons with momenta fci, k 2 — from the 
external legs: 



M^(p,q) 



— i 



i — 



(24) 



d v G~ x 
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Inserting (|24| ) into (^) and replacing g(0) by g{k) we get 



d 2 T»(p,q) = 

= g{q){A v {q 2 )d v Y»{p, q) + d u V*(p, q)A v ( qi ) 
A u (q 2 )T^p,q)A v ( qi )} + d 2 T^p,q). 



(25) 



Here we have introced the notations 



ft l9 = ?±|, = d ll G~ l {q)G{q) and Ap(q) = G(q)d lt <r 1 (q). (26) 



The correction terms c? 2 f M are defined as a set of diagrams with exact Green's 
functions and vertices. They contain first order derivatives of both the photon 
lines and the positron lines, second order derivatives of the positron lines and, 
in the same way as in the case of the Green's function, corrections due to the 
replacement of g(0) by g(q). 

The same equation is valid for the interaction amplitude of fermions with the 
external field provided this interaction does not depend on the relative momentum 
q of the fermions. 

The equations for the interaction with external fields are essential. Indeed, if 
these equations have solutions that decrease at large virtualities of the fermions 
- i.e. solutions which do not require driving terms - this means the existence of 
bound states. 

For the sake of convenience we rewrite the equation fl2"5p in the form 



d 2 <j)(p, q) = g(q){A v (q 2 )d v (f>(p, q) + d v <j>(p, q)A v ( qi ) - A v {q 2 )(j){p, q)A I/ (q 1 )} (27) 



which will be understood as the equation for the bound state with a spin which 
is defined by the invariant structure of the matrix <fi. It is important to note 
that the accuracy of the equations for the vertex T^(p, q) and for the bound 
states differs from the accuracy of the equation for the Green's function. The 
functions T^(p, q) and 0(p, q) depend, among others, on the ratio p A jq\q\- If this 
parameter becomes large, then Y ^ and 4>(p, q) contain, in general, the so-called 
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Sudakov logarithms a In which are not included in the equations (Pq) and 
(j2"7D . Hence, the equations are valid only if 



4 Equations for Green's functions in QCD 

QCD is the theory of interacting quarks and gluons. The description of quarks is 
more or less the same as in QED. Gluons, however, are very different. Even the 
fact that being a spin 1 particle the gluon has to have a multi-component wave 




(28) 
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function is not seen explicitly. In the usual approach it is described by a Green's 
function D^ u = (A^x) , A v (y)) which in momentum space (in Feynman gauge) 
always can be written in the form 

= ^C(k 2 ). (29) 

It contains only one unknown function. All spin properties of the gluon are in- 
cluded in the momentum dependence of its interaction. In order to introduce 
multi-component Green's functions, we have to formulate the theory in a form in 
which the interaction is momentum independent: we have to replace the usual 
description of the gluons and their interactions by a Duffin-Kemmer type formu- 
lation. In the framework of this description the gluon Lagrangian is 

C(x) = {d„A v - dvAp + g[A u , A^F^ + (30) 

The potential A^ and the field strength F^ u are independent quantities. The 
interaction is momentum independent and equals [A^, A^F^. In this formulation 
we have three independent Green's functions 

(A^AJ, (F^,A P ), (F^F^). (31) 

In order to fix the gauge in a covariant way, it is, of course, necessary to introduce 
ghosts by adding a gauge-fixing term 

AC= C -(d,A,) 2 . (32) 

The three Green's functions can be combined into one by introducing the ten- 
component state 

We use the parameter m, which has the dimension of mass, to convert the lower 
component of \I/ into the same dimension as the upper component. The equations 
for the fields corresponding to the Lagrangian (ftOD, (j3"2"D are 



+ C^d v A v = 0, 

d lt A v -d u Ar + g[A lt ,A u )-F lu/ = 0. (34) 
In terms of the state \I/ we have 

/3 M H<9 M + gA^ - m 7 _^ - ^-Q{p 2 - p 2 )^ = 0; (35) 
V 2 = -PpPudpd u , p 2 = <9 M <9 M . 
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Pn are Duffin-Kemmer matrices satisfying the commutation relation 

A Ac A + A Ac A = Sikft + Supn 



(36) 



they are connecting the upper and lower components of \I/ and have a 
simple representation 



(/3 



o"7 



l 

71 



( ^/icr *5p7 ^7 ^pcr ) • 



(37) 



The quantities 7± are projector operators for the upper and lower components of 

If we want to preserve the dimension of the (A^, A u ) component of the free 
Green's function, D has to satisfy the equation 



m m 



1 

m 2 



The solution of this equation is 



(38) 



(39) 



In Feynman gauge £ = 1 and C\ — 1, in Landau gauge = The three terms 
in (|3"5|) correspond to the three independent Green's functions (|3T|) . 

The vertex for the interaction of three gluons with momenta Pi,P2,P3, colours 
a, b, c and Duffin-Kemmer indices a, A 7 has the structure 



a,a,pi;6,/3,p2;c,7,p3 



P2 



P2 



Pi 



+ 



P3 



Pi 



+ 



P3 



Pi 




P2 
P3 



(40) 



/' 



+ 



/' 



P" 



v pa pa v 

The coupling constant remains in our notation g. 

Let us consider in this approach the properties of the exact Green's function 



IT 1 = k - 7- - ((k 2 - k 2 ) - S - <5C(A: 2 - fc 2 ^ 



(41) 



P 



where £ is defined digrammatically: 





(42) 
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It contains four matrix elements E ++ , E |_, E H , E It can be also represented 

in the form 

E = fcE 1 +7_E 2 + 7+ A; 2 E3. (43) 

The factor k 2 in the third term is necessary to preserve the current conservation; 
in first order, the ghost contributes only to E 3 . 

We added the term 5( in (^Tj) to be able to fix the gauge for the exact Green's 
function. Instead of (Ell) we can write 



(44) 



zr 1 = z~ l k - Z'S- - z~ l k 2 - c(k 2 - k 2 ). 

The Green's function then will be equal 



D 



(ZT* - Z 2 'Z^)p 2 



where 



£7 = 1 + 



z 2 z x 



- 1 



In Feynman gauge 



C — Z 2 Z X z — Z 3 X . 



(45) 



(46) 



(47) 



In order to understand the meaning of the denominator in (]4"5|), let us consider 
the renormalization properties of simple diagrams, for example E + __ + E |_. 




r ++ - + (+->-) 



(4£ 



where g^m is the effective bare coupling; 

r+-+ ~ zf . 

If this is true, we have in Feynman gauge 



Z^Z 2 X 



{z^-z 2 i z^r 



(49) 



(50) 



However, p M <9 M E has to be proportional to aZ l 1 . This means that we have to 

expect 

(51) 



z^z 2 l 

f ry — 2 ry — 1 ry- 
v^l — ^2 ^3 



a (Pl z -l 



Q'O 



i.e. 



Z 2 Z^ 2 — Z 3 1 — */ — Z\ 2 Z 2 . 



a 
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This is our definition of az(p). It has all known properties of the renormalized 
coupling a. As a result, D can be written in the form 



D = J^-ZlZl{Z 2 Z^k + C 1+ + (Z 2 Z^k) 2 ^}\ + ^t- (52) 
V ocq p m z 

In this approach ao is not a quantity coming from the normalization: it is the 
bare coupling. The theory has to be defined as the limit a® — > 0. In this context 
the expression (|52| ) has a very interesting property. In the limit «o — > the 
Green's function D contains only Z^ 1 and Z 2 l . According to (|5T|) , in this limit 
Z^ 1 = Zi 2 Z 2 . Because of this, when we regard the equations for Zf 1 and Z 2 X in 
the way we did for the fermionic Green's function in QED, ao has to disappear. 
Consequently, we have equations for Zf 1 and Z 2 X with a(p 2 ) being arbitrary. 
The equation for Z^ 1 will not help since, due to the equality Z 3 = Z^ 1 Z 2) it 
has to be an identity. We will see that an equation for a appears when we will 
consider the correction of the order of y/ao. 

Before formulating the equation for the Green's function, let us see what 
the Ward identity looks like in this formulation. Consider the relation between 
(p 2 — Pi) l i^ tJ '{p2,Pi) and the Green's function for the bare vertex V 1 = f/3 11 : 

P^Pn =fo~Pi = rn[k 2 - 7_ - (fa - 7-)] = (53) 

= m[D-\k 2 ) - D-\fa)} + m ( [(k 2 2 - k 2 2 ) - (k 2 - k 2 )} (54) 

At p -> 

rj| p=0 = [md.D- 1 + m( dXk 2 ~ k 2 )}f (55) 

which is, of course, the usual complication due to the Slavnov- Taylor Ward iden- 
tity. But in this formalism the additional term is equal 

Introducing this vertex in an arbitrary diagram we obtain the relation 

r M = m^ZT 1 - & JLfyD- 1 ]. (56) 

According to (fill), (|47"|) , in Feynman gauge we can write 

r M = /m^[(zr l -a c zr i )^-7_(z 2 - i -a c z 2 - i )-7 + ^ 2 (^r%-a c (zr 2 z 2 ))]. (57) 

The quantities c^Z-f , d^Z 2 x must be calculated from the equations for Z^ 1 
and Z 2 l . But if the equations are formulated in terms of the exact Green's 
functions, the dependence on ( enters only through these Green's functions D 
which, according to (f45|), (f46|), include ( only through the quantity C. In the 
limit «o — > 0, however, C does not depend on (. Consequently, d^Z{ 1 and d^Z^ 1 
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are equal to zero and we have the following simple relation for the vertex at zero 
momentum: 

T li = fmd ll b- 1 . (58) 
D^ 1 is defined by (0); it does not contain gauge-fixing terms. 

L)- 1 = Z^k - Z 2 " 1 7_ - 1+ Z^ 2 Z 2 k 2 . (59) 

Let us first consider the equations for Zf 1 and Z^ 1 in the same way as we did 
for fermions. As in the case of QED, we will use the Feynman gauge in order to 
simplify the one-gluon contribution to the equation for the Green's function. 




P~Pi P-P2 . 

+ KS> (60) 

Pi P2 

Each line here contains the exact propagator flBlq). Differentiating the propagator 
along the upper line, we obtain the contribution corresponding to the second 
derivative of one line 

d 2 D = -4n 2 i5 4 (p) 1+ Z^0) - ^fd.NZ^ + \d 2 NZ^ + d 2 ^- (61) 

jr p z m 2 

where 

Z = Z 2 Z 2 ^ (62) 
a 

and we introduce the notation 

Zz Zol- 

D = —N+-^. 63 

The contribution of this derivative to the equation will have the form 

- 4n 2 iZ^M(p,p) + M l . (64) 

Here M(p,p) is the gluon-gluon scattering amplitude at zero momentum of one 
of the gluons. The second term Mi is defined diagrammatically: 

M 1 = 



Taking the first derivative of two different lines, we get the contribution 

M 2 = -<^^>- + -^J^>- + .... (65) 

As a result we have 

d 2 (Z{ 1 k - Z^y-) = —Z^M{p,p) + Mi + M 2 . (66) 

7T 
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We have to remember that on the right-hand side we have to take only the matrix 
elements ( — (+|— ) and (— |— ). Writing 



M{p, p) = / 2 r M Dr" + M(p, p) (67) 



and using the Ward identity fl58|), we obtain an equation of the same structure 
as for fermions: 

= f^l 2 ZH0)zHp)d,D-^N + ^^]d^ + L-, (61 

I. M ■ M. ■ A I , , f 2 = n c = 3. 

If we now make the same trick as before: replace Z(0) by Z(p) and redefine L, 
we will have in the limit «o — > 

d 2 Z^, k] = 3^1[ 7 _, d,D-\z{k + Z\Z^ 1+ + Z.k^d,!)- 1 ] + [ T _, L'], 

7T K 

(69) 

c^V = Z^L±- 1 _d^D-\Z l k + Z 2 Z^ 1+ + Z 2 fe 2 7 _)^^ _1 7- + 7-^7- 

(70) 

The equation for the fermionic Green's function in QCD will differ from (^) only 
by the factor A a A a = | (A a are colour matrices). The reason is the following. The 
derivation of the equation fl2"T| ) was based on the relation between the fermionic 
Green's function and the amplitude of zero momentum photon emission by a 
fermion 

T ll (q,0)=d ll G- 1 (q). (71) 

In the usual formulation this relation is not correct in QCD. The simple relation 
( |58|) for the amplitude of the zero momentum gluon emission by a gluon implies, 
however, that the amplitude of a zero momentum gluon emission by a quark has 
to be equal 

T«( q ,0) = \ a (Z^Z 2 ) l *d,G- l (q). (72) 

Together with (|5^ ) it leads to the equation (^Tj). The equation for the colourless 
vertices remains also the same. The equation for a three-gluon vertex, however, 
will be essentially different. In the same way as for the fermionic case, but taking 
into account the non-commutativity of the gluon coupling, we can show that 

r = ir^) = f M 

x upper lA 1 jap 1 

satisfies the following equation: 
d 2 t» (q + Z,p, q -E 
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A v = d v D- 1 D , ~A v = Dd v D- x - p x = ? + | , P3 = g-|- 

The right-hand side corresponds to all possible gluon emissions from the external 
lines: 

1 d„b- 1 



+ 



d v D- 



+ 



d u D~ 



+ 



+ 



d v D- 



d v b- x r d v D~ x dvb- 1 dub- 1 r 

Higher order terms can be written in the same diagrammatic way. The third 
term in (ff^) and in the corresponding diagrammatic expression equals zero at 
a = since DTD = at a = 0. 

The equation (|73|) is, indeed, quite different from the equation for a colourless 
vertex. The main difference comes from the fact that on the right-hand side it 
contains derivatives of T not only over q but also over p. Hence, in order to find 
T, we have to write three different equations for second derivatives over three 
different momenta. 

A similar equation is valid for the quark-gluon vertex. 

Knowing the Green's function and the vertices, all the other amplitudes for 
the interactions and quarks and gluons can be written in the usual perturbative 
way. These amplitudes have no divergences and contain inside the gluon Green's 
function the unknown function a(p). To formulate the theory in an unambiguous 
way, without any references to the cutoff and the regularization, we have to find 
the equation for a(p) as we did in QED, and learn, how to write the higher terms 
more elegantly and constructively. I postpone the investigation of this problem 
to another paper. What I want to discuss now is the main difference between an 
asymptotically free theory and an infrared free theory. 

The equations for Green's functions of quarks and gluons are proven to be 
second order integro-differential equations. To solve them, we need boundary 
conditions. In an infrared theory the boundary conditions are known: they are 
the conditions for the existence of free elementary particles at small momenta. 
In an asymptotically free theory the interaction is small at large momenta, and 
we expect to have here a perturbative solution. However, in the region of large 
momenta all the equations have two types of solutions: the perturbative solution 
and solutions which decrease as some power of the momenta and therefore contain 
dimensional parameters reflecting the density of different condensates. These 
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parameters have to be defined either by the introduction of additional conditions 
of the type of conservation laws or by the self-consistency of the solutions in the 
small momentum region (or both). The next section of this paper will be devoted, 
essentially, to the discussion of this problem. 



5 Spontaneous symmetry breaking in asymp- 
totically free theories 

In the present section we will consider the equation for the fermion Green's func- 
tion in QCD 

d 2 G-\q) = g(q)d,G-\q)G(q)d,G-\q) (74) 



where 



g(q) = °m A» A « 



7T 



Aa(q 2 
3tt 



This equation was discussed extensively [1] in connection with the problem of 
quark confinement. Here we shall use the equation for the discussion of the 
spontaneous breaking of chiral symmetry in asymptotically free theories. The 
asymptotic freedom is reflected in this equation by the fact that a(q 2 ) decreases 
when q 2 — > oo. We will assume that in the limit q — > a(q 2 ) approaches a finite 
value. 

9 




Figure 1 



At q — > oo the solution of the equation (|74j) has the form 



G-\q) 



[m ~ q) + ^ + "T 



(75) 



q = IiiQim- 

If a = 0, the quantities Z, m, v-i are arbitrary parameters. In the case when 
a(q 2 ) is defined by perturbation theory, Z,m, z/ x , u 2 are the following functions: 



Z = Z 



a \ 
a J 



a\~< m n fay 1 - 2 

m = m [ — , v x>2 = ^12 — 



(76) 
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The anomalous dimensions 7^, 7 m and 71,2 can easily be found from the equation 
(|74|). Generally speaking, the solution depends on four parameters. In the limit 
q 2 — > 00 the chiral invariant solution can be written as 

G- 1 = -Zq (l-fj. (77) 

The general solution ([75|) corresponds to massive quarks. In the solution which 
corresponds to spontaneously broken chiral symmetry, m Q = 0. In this case the 
mass term decreases when q 2 — > 00; the term v 2 is responsible for the violation 
of the symmetry. 

Multiplying (74) by G(q)G~ 1 (q) we obtain the equation 

d 2 G-\q) = gA,{q)A,{q)G-\q) (78) 

where 

A^q)=d,G-\q)G{q). (79) 

Clearly, it has a structure which corresponds to the equation for particle propa- 
gation in the self-consistent field gA^A^. It is easily seen that ( |74"D is equivalent 
to the equation for A^ of the form 

d^ = -/LW (3=l-g. (80) 

The matrix G' 1 is defined by two invariant functions and can be written as 

G- 1 = Z-\q)[m(q) - q) = pe~^ (81) 



where 



A^ is here 



- q 

n = -. 

q 



Inserting fl8~l|) and fl82| ) into (|78|) or (|79"D, we obtain a set of non-linear equations 
for p and 0. We can linearize the equation for p by writing 

-(f)* «) 

for a constant /3 or 

p = _ e JiU- 1 J7 (84) 
for (3 which is a function of q. Here 

u = q u d u u = ^- } £ = lng. 
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As a result, we get for u and the following set of equations: 



u-u + (3 2 



3 sinh - + — 
2 4 

3 sinh . 



u 



+ 2-0 -3 sinh = 0. 

u 



(85) 
(86) 



For a constant (3 the conservation law 

d^E = 

is fulfilled; 



E = u 2 -u 2 + {3 2 



3 sinh 



u 2 . 



(87) 



The term ^{u — u) is of the order of g 2 and it can almost always be neglected. 

The asymptotic behaviour (|75|) of the Green's function in the limit q 2 — » oo 
(/3 — > 1) corresponds to 

u — > Cq 2 , — > in. 

The chiral invariant solution corresponds to = zV. Close to the value 
(i.e. = i7r + 0), at large g 2 we have 



/7T 



1 + 3/? 2 . 



Hence, ( ^6|) is an oscillator equation with damping if q 2 is growing and with 
acceleration if q 2 is decreasing: 



+ 



2Jl + 3^ 



19) 



This means that the chiral invariant solution = in is unstable in an asymptot- 
ically free theory. 

Let us consider the equation for at negative q 2 values in detail. Due to 
fl8i"D, = iif} is in this case purely negative and ( p6|) describes the motion of the 
particle as a function of the "time" £ in a periodic field; the damping (or the 
acceleration) is defined by (|88|). 
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At £ — > oo (g 2 — > oo) the particle is situated at one of the minima of the potential; 
it accelerates as £ decreases. The acceleration rate is defined by the parameters 
of the solution (p3|) in the limit q 2 — > oo. If q is decreasing we have, generally 
speaking, two possible behaviours for the solution. It goes to infinity if - remains 
positive all the time, or it may approach again a minimum of the potential if - 
changes sign. In the latter case it is easy to see that G^ 1 has a singularity as 
q -> 0. 

There is only one possibility to avoid having a singularity in the Euclidean 
region including q — 0. We have to choose the parameters which determine the 
acceleration at q 2 — > oo so that the particle appears at the maximum of the 
potential if q — > (£ — > — oo). In order to find exceptional solutions without 
singularities at q 2 < it is natural to solve the equation by fixing the solution at 
q — > (£ — > — oo). The solution of (jS5[),(jS6D corresponding to a maximum (e.g. 
■0 = 0) at g — >■ is 



iip = — , M = Z" 1 (0)g (90) 

rrir 



m 

In the case of such a solution the constant E equals zero, and 

(91) 



u 
u 



\ 



l + P 2 



q • 2^ . ^ 

3 sin — H — - 



Inserting fl9T|) into (p6|) we obtain one non-linear equation. It can be analysed 
easily for arbitrary q 2 values. The solution which we are interested in contains 
essentially one parameter m c which can be related to the renormalized fermion 
mass. The parameter Z~ x is irrelevant; it defines the normalization of the Green's 
function at q = and can be chosen as 1. The solution of the equation leads 
to the unambiguous determination of the asymptotic parameters of the Green's 
function ( |75|) by m c and by the parameter A of the strong interaction which enter 

As it was mentioned before, to the spontaneous breaking of chiral symmetry 
corresponds an asymptotic behaviour of G' 1 in which m — 0. The existence 
of such a solution requires a connection between m c and A. The renormalized 
fermion mass as well as the other "condensate" parameters are then determined 
by the strong interaction parameter A. 

Let us consider this solution in detail. Starting from the point ip = at 
£ — > — oo, the solution will, obviously, reach the minimum of the potential either 
monotonically (if the damping is strong enough) or in an oscillating way. In our 
case the damping depends on the value of (5. If g is small, (3 is close to unity, the 
damping is strong and the solution has a monotonic behaviour. By decreasing 
P the solution may become an oscillating one. For obtaining the value of (3 at 
which the oscillation begins there is no need to solve the equation at any q. It 
will be sufficient to investigate the solution in the region where ip becomes close 
to 7r; here ip = ir + ip and ip satisfies the equation (B3). The solution can be 
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written in the form 



ij, = e- p? Ccos(^2 - 3f3 2 £ + 5), 
which oscillates if 



V 



1 + 3/? 2 



^<3 5 1 -V3 < ^ <1 + V3- 
Oscillations in ip mean that the mass term in (j8l|) 

i \ 

m(q) ~ i— q 



(92) 
(93) 

(94) 



starts to oscillate. By solving the equation for bound states we can check that 
at g > g c bound states appear with wave functions behaving like fl92]) which 
coincides with the behaviour of the solution of the Dirac equation in the field of a 
point-like static charge Ze when Z > 137. The simplest example for such bound 
states are Goldstone states the wave functions of which, as we shall see it in the 
next section, are proportional to m(q). 

We have found the oscillations and the critical value g c — 1 — using the 
assumption that g = const. In reality g depends on q (see Fig.l) and the situation 
is somewhat more complicated. It reminds the case of the equation for a critical 
charge of finite radius. In the region of small q values where g(q) is close to a 
constant g(q) ~ g(0) we can consider if) as an independent variable and q 2 as a 
function of ip. It can be shown that for g(0) satisfying the condition (|95|) there 
are two regions < ip < tt — ip c , n + ip c < ip < 2n in the q 2 , ip plane (Fig.3) 




Figure 3 

where the solution if)(q) is a monotonic one and there is a region 7r — ip c < ip < 
n + tp c where the solution oscillates. The value of ip c is determined by the equality 

sin 2 f = ( = - 2 ) \r~r^T- 7= ~— = ; (95) 
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it differs from zero if f3 2 < |. Considering f3 as a function of g 2 in fl9"5| ) and taking 
into account that (3 2 — > 1 at q 2 — > oo, we obtain a region inside the dotted curve 
where the solution oscillates. There are no oscillations if q 2 > X 2 (j3 2 (\ 2 ) = |). 
Due to (ff5|) and (|76|), we can write in the region q 2 ^> A 2 



J(^-t) = — (-)' + ^(-)'"- 06) 



2 g \a / g 3 V«o 

If /? 2 > |, the solution which equals ip = at g — > transforms monotonically 

into (j96f) where m 7^ 0. If /5 2 < |, ip{\ 2 ) and V'(^ 2 ) start to oscillate as functions 
of m c , and mo can turn into zero at a certain m c value. This means that we have 
a solution corresponding to broken chiral symmetry. 

If mo = 0, there exist also a large number of solutions depending on the 
parameters Ui, v?,- With the same sign of V\ we can have different solutions 
corresponding to the curves I and II which in the limit q — > reach if) = and 
ip = 2ir, respectively. The solutions ^(0) = 2tt correspond to smaller values u±, 
m c . 

Let us consider the solutions of types I and II in detail for complex q - this 
is justified since G _1 has to satisfy the requirements of analiticity and unitarity. 
We will show that both solutions have singularities at real positive g 2 values. The 
solutions are chosen in such a way that they are regular as q — ► and have no 
singularities at g 2 < 0. Due to the analiticity of the equations, the behaviour of 
the solution for g 2 > can be found by solving the same equations (|8~6p, fl91)) 
with the same boundary conditions at q —>■ 0. 



If (3 is fixed, the equations (0), (0) can be rewritten in a simpler form. 
Denote 



- =p(<f>)- 
u 



Then 



|^ = -/W + 3sinh 2 ^-l (97) 



^ = \v +3sinh2 2 _1 (98) 

with the boundary condition p = 1 at = for a I-type solution. The phase 
diagram corresponding to the equation (|97D is shown in Fig. 4 where the solid line 



represents the solution of the equation 

p 2 = l-3sinh 2 -. (99) 
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The function p = p(4>) has the shape of the dashed line in Fig.4. 

kp 




Figure 4 

The ^-dependence of p at —* oo is different for j3 > \ and (3 < |. In both cases 
approaches infinity at finite £ values. We have at (3 > ~, — ► oo 



C > 



(100) 



and at /3 < ~, — > oo 



p = -2/3 2 C ie ^ 



2 V 1 - 4/? 2 

The equation fl98|) enables us to find £ as a function of 0. 



(101) 



p 2 (0) + 3sinh 2 f - 1 



(102) 



The integral in (|102| ) converges in both cases {(3 > | , /3 < |); because of this, 
goes to infinity at a finite £ = £*(<? — * m*). Near q = m* at (3 > | (g < |) we 
have 



u = tio\/l 

m 



e 2 



(103) 



and, consequently, 



q \p 1 ,^ , / 1 \ p 1 



G^fe) = Z?\ ( 1 - JL) ' + g) + y -(<?-<?)}. (104) 



If/?<§ (g>§), 



W = ( 1 — 



Ci 1 



_9_ 



(105) 
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and thus 



q \ 2 1 1 



(106) 



It can be easily shown that the relation between the position of the singularity 
q = m* of the Green's function and the quantity m c can be written as 



In 



m 



P 



2\/p 2 (0) + 3sinh 2 | - 1 sinh20 



(107) 



The formulae (|104|) and (|106|) define the behaviour of G~ x (q) near the singularity 
for a solution of type I. In order to obtain the behaviour of G _1 near the sin- 
gularity for a solution of type II it is sufficient to notice that at q — > such a 
solution has the form 

= 2m - — . (108) 



m' 



The replacement of by + 2iri changes only the sign of G 1 . Replacing 
by — and solving (P71) and (R3) for q > we will find to be negative, and 
near the singularity (f) will go to — oo. This means the following. If the first term 
G+^g) °f the expression 



G-\q)=G-\q)- 1 



(109) 



equals zero for the solution I at q — m*, then G__ l is zero for the solution II at 
q = m*. If q 2 > m* 2 , both solutions become complex: 



Z7T 1 

+ 7^ ' 0> 2 



-2 hid 



Q 



1 + 2ivr 



/3< 



1 



(110) 



The trajectories of 0(g) at < q < oo in the complex plane are shown in Fig.5 
for the solutions I, II. 



II 



ir 



2Z7T 



■ITT 



Figure 5 
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The solution of the type I has a remarkable feature: 30 > it for any q > m* 
values. Hence, taking = 0i + ifa, the imaginary part of m(q) in (|81[) 



g CV COSh 2 ( 

sinh | ( 



+ ic 



-g sm 02 



2|sinh§( 



+ H 



(111) 



turns out to be positive. At the same time 5sm(q) for the solution of the type 
II is an oscillating function; this can lead to a contradiction with the unitarity 
condition for the Green's function. For complex q values G~ 1 (q) has no singular- 
ities. Moving in the complex plane along an arbitrary ray the trajectory of 0(g) 
doesn't approach infinity (curves I' and II'). 

The main result of this section is the following. In the framework of the 
equation for the fermion Green's function (|74j ) there exist solutions corresponding 
to broken symmetry provided g(q) has an asymptotically free behaviour and 
g(0) > 1 — •*/-. These solutions behave at q 2 — > oo as 



q 1 



(112) 



The expression (|112) has a mass term decreasing at infinity. 



6 Axial current conservation and Goldstone 
states 

If a fermion Green's function corresponds to symmetry breaking, it is natural 
to expect the existence of Goldstone-type bound states. This expectation is 
connected with the belief that if mo = 0, the axial current has to be conserved. 
This is, however, not necessarily true because due to divergences in the theory a 
leakage of the current is possible in the region of the ultraviolet cutoff. A typical 
example for this phenomenon is the anomaly. But even in a non-anomalous case 
it is not obvious whether the current conservation is contained by the equation 
for the Green's function or it is imposed as a condition on the solution of the 
equation. In order to clarify this, let us consider the equation for the bound state 
supposing that it is a pseudoscalar. 

<9 2 0(p,g) = g(q){AMd^(p,q) + d^,q)A^ qi ) - A lx {q 2 )(j ) {p,q)A ll {q x )} (113) 

This equation has to have a solution decreasing at large q 2 and p 2 = 0. It is easy 
to see that indeed there exists such a solution. At p = the equation for is 
an equation for the variation of a function which satisfies the equation for G^ 1 . 
If this variation is taken in the form = C{75,G -1 }, obviously satisfies the 
equation and decreases as q 2 — > oo 

-> C^ 7s . (114) 
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This, however, does not mean that we have particles with p 2 = 0. Indeed, the 
equation (\113\) has a solution decreasing at q 2 — > oo for any p 2 values. The reason 



for this is that the equation is highly degenerate. It has a solution of the form 

= O 1 G-\q 1 )+G-\q 2 )O 2 (115) 
where 0\ and 2 are any combinations of Dirac matrices with coefficients not 



depending on q. This can be checked by substituting directly ( |115| ) into 
and using the equation for G _1 .The reason for this degeneracy is the invariance of 
the equation ( |1 13| ) under Lorentz transformations, under all the possible discrete 



symmetries of the Dirac equation and for a fixed g = const even scale invariance 
and translational invariance in momentum space. If, for example, 0\ = 2 = 75, 
then, due to the Ward identity 



p,Tl(p,q) = l5 G-\q 1 ) + G- 1 (q 2 h 5 , (116) 

<P = l 5 G- 1 (q 1 )+G~ 1 (q 2 ) l5 (117) 
is the divergence of the axial current. Hence, if the equation ( |1 13|) is satisfied for 



PfiT 5 ^ it has to have the solution ( |116| ). 

In order to show that ( |113j ) has a decreasing solution at q 2 — > 00 for any p 
value, let us notice that in Euclidean space this equation has the structure of the 
Schrodinger equation 

{-d 2 + U)^> = e^> (118) 

in four dimensions at e — 0, with a potential depending on q, spin variables and 
the external vector p. For such a potential the total four-dimensional angular 
momentum is not conserved, only its projection fi onto p v . An equation of this 
type has always non-singular solutions with incoming waves of given [i. E.g., 
considering an incoming wave with fi = for the pseudoscalar <ft, we will have a 
solution 

00 = I5C0 + decreasing scattered waves 

at q 2 — > 00, or, regarding an incoming wave of the form 75^7^ = 75P, we will 
have 

01 — I5P + decreasing scattered waves. 

Suppose we found the solution <pi. In this case, due to the fact that for q — > 
the solution (|116|) behaves as 

l5 G- 1 (q 1 ) + G-\q 2 ) l5 - Z ( l5 p + ^ ) , (119) 



Q 2 



we will find that 



4> = Z<t> x - l 5 G-\ qi ) - G-\q 2 ) lb - (120) 



q 2 
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is decreasing with q 2 — > oo at any p. 

By stating that the equation for bound states has a solution at any p values 
does not mean that there are no bound states; it means only that the mass of 
the bound state has to be calculated independently The most natural way to do 
this is to calculate the self-energy of the state 



E(p) 



and to solve the equation 




E(p') = 0, 

or to calculate the forward Compton scattering 



(121) 



(122) 






(123) 



of the bound state on a fermion and after that integrate over the distribution 
of fermions in the vacuum. But this way we will never get a massless Gold- 
stone state. The reason for this is the almost obvious fact that ( |121| ) and ( |123j ) 
contradict the condition of axial current conservation. 

Let us consider the condition for current conservation in detail. If we introduce 

as a set of diagrams 



7^75 + 



7^75 



+ 



(124) 



with a massive fermionic Green's function, it will not satisfy the Ward identity. 
However, the Ward identity will be satisfied by the sum of and of the Goldstone 
contribution 




+ 



Q2 qi 




(125) 



l 5 G-\ qi ) + G'\q 2 ) l5 . 



(126) 



Here g is the Yukawa coupling of a Goldstone boson to a fermion: 



ifPiA = 7^75 V^^/ 9 ' 
27 



(127) 



Knowing p^T 5 ^, we can define the Yukawa coupling g by ( |126| ), ( |127| ). Let us 
apply the operator — d 2 + U to (|126|) . We find that gf satisfies the equation ( |113| ) 
since p M r^ and the right-hand-side of ( |126|) satisfy the same equation. But the 
existence of ( |125|) implies that the mass of the Goldstone boson has to be zero. 

In order to clarify the situation with Compton scattering, let us consider the 
Ward identity for the amplitude 



rl(k',q 2 , qi ,k) 




(128) 



where the dotted line with the momentum k' corresponds to the axial current, 
the wavy line with k corresponds to the Goldstone state and the other two lines 
to the fermions. The Ward identity for this amplitude is 



k' 



k'^l(q 2 , qi ,k) 



75 



g(k,qi) 





life' 



(129) 



9(Q2,k) 



75 



To fulfil ( |129| ), T^(g 2 , qi, k) has to be equal 



rj(g2,?i,fc) 



k! 
r 5 



k' 



g(qi,ky 
72 qi 



+ 



r 5 

A 1 , 



Q2 qi 



+ 




; (130) 

A(<?2<?1, k', k) 



here is defined by (|125|) . Using ( |125|) and ( |129| ) we obtain 

«/A = 759(k, qi) + g(q 2 , k)^. 



At large qf, q\ values we will have 

.275^ 
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A 



4z/ 3 

7v 



(131) 



(132) 



This means that similarly to the case of asymptotically non-free theories the 
Goldstone - fermion scattering amplitude does not depend on the momentum of 
the Goldstone boson; it decreases only as a function of fermion virtuality. Under 
these circumstances it is obvious that even in an asymptotically free theory the 
Goldstone boson has a point-like structure. 
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Amplitudes for the interactions of many Goldstones with fermions can be 
found in an analogous way and have the same properties. The self-energy of the 
Goldstone state is now different from (|121| ). It contains two terms 

m = — + jQ^ (133) 

and equals p 2 at small p 2 . 

As we already have said, in this approach the Ward identity becomes the 
definition of the Yukawa coupling g (the wave function of the Goldstone boson) 
through p^T 5 ^ which has a clear diagrammatic meaning. The equation contains 
the amplitude / of the Goldstone - axial current transition. The expression fll27| ) 
for fpp is highly symbolic, because it contains overlapping divergences. In order 
to write a sensible expression we can use the same procedure as we did in section 
|2| when we calculated the polarization operator of the photon. Applying the Ward 
identity, we can write instead of ( |127p 

f% = 7m75<^^> I>„ • (134) 
Differentiating (|3|) over p we will get for p = 




(135) 



If we want to get rid of the 75-s we have to commute 75 with the Green's functions 
and the interaction vertices along one of the fermionic lines. Doing this, we obtain 

{75,G" 1 }75 {75, G- 1 } { 75; G- 1 } 




(136) 



Due to the conservation of the vector current, the first two terms in (|136|) are 
zero at p = 0. The first one is just the photon polarization operator at p = 0, the 
second one is the amplitude for the decay of a zero momentum scalar into two zero 
momentum photons. The last term looks like the zero momentum pseudoscalar 
- photon scattering amplitude which also has to be zero. This, however, is not 
true, because it does not contain all the necessary diagrams. It does not involve 
overlapping divergences. As a result, we can write 



{75, G- 1 } {^G- 1 } 



{75, G- 1 } ^G- 1 } 



4/ 2 



d.G" 1 + d,G- 1 



d u G- L + 
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^G- 1 } {75, G- 1 } 
+ d^G- 1 <"^> d^G- 1 + ■■■ . (137) 



In the zeroth order of - 

7T 

f = \S (^ Tr{75 ' G ^ G ^ G-^GA^Mq)- (138) 

7 Flavour singlet and flavour non-singlet Gold- 
stones states. The U(l) problem 

Up to now we have discussed the Goldstone states in a relatively abstract way 
without fixing the concrete asymptotically free theory. In real QCD we have 
quarks with different flavours and there is a difference between flavour singlet 
and flavour non-singlet states. In order to clarify the picture it will be useful to 
describe the Goldstone state in a different way. 

The previous discussion shows that the Goldstone states in asymptotically 
free and non-free theories are rather similar. Therefore it is natural to try to 
introduce the Goldstone boson in the usual way as a point-like state and to see 
how this state will interact. In the usual discussion of a Goldstone particle we 
suppose that there is a point-like pseudoscalar interaction between this particle 
and a fermion with the pseudovector coupling 

Sp i 

< = P^ T - (139) 
« \ « Jo 

This interaction induces the radiative correction to the propagator D{p 2 ) of the 
pseudoscalar 

D{p 2 ) = vww + 'vK^^^yv^ H ; (140) 

Dq is the bare pseudoscalar Green's function. If the fermion is massless and the 
axial current is conserved, this pseudoscalar will not interact; its self-energy 

nv 2 ) = 4^7^75^ (141) 
Jo 

is equal to zero. If, due to symmetry breaking, the fermion becomes massive, 
it starts to interact and acquires a self-energy different from zero. By using the 
diagrammatic definition of we will find 



= ^ (142) 
Jo 
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where / is the same amplitude for the Goldstone - current transition as what 
was discussed in the previous section. Hence, 



D( P 2 



f 2 
Jo 



2 f2 



In the limit fo — > we will have 



D(p 2 



/o 2 



(143) 



(144) 



the pseudovector Goldstone - fermion interaction is jPpX^ with a pseudovector 
coupling j defined by the fermion mass. The limiting procedure fo — > can be 
understood if we accept that the interaction responsible for symmetry breaking 
changes the fermion vacuum fluctuations not only at finite momenta but also 
near the ultraviolet cutoff. This change in the fermion vacuum fluctuations is 
responsible for the leakage of the axial current in the region of finite momenta; 
it can produce the driving term for the Goldstone state, recovering the current 
conservation. 

In general, the pseudovector coupling has a disadvantage compared to the 
pseudoscalar coupling which we have discussed before: it looks unrenormalizable. 
But in the case of flavour non-singlet states it can always be replaced by the 
pseudoscalar coupling with the help of the trivial relation 



P75 = (<?2 - m (<?2))75 + 75(gi - rn(qi)) + m(g 2 )7s + 75^(?i 



(145) 



which leads to the Ward identity ( |126| ) for pseudoscalar coupling. If we include 
the emission and the absorption of Goldstone bosons inside the diagram, then in 
the process of this replacement a point-like amplitude appears which corresponds 
to the quark interaction with many Goldstone bosons. Nevertheless it is possible 
to prove that, due to the decrease of this amplitude as the function of quark 
virtuality, the theory is renormalizable. 

Due to the anomaly, the case of a flavour singlet current is very different even 
on the level of the Goldstone Green's function. In this case the corresponding 
polarization operator E will contain not only the quark loop which we have 
discussed but also a gluonic contribution 





(146) 



The triangle diagram f^ v included in ( |146| ) was calculated many years ago by 
Adler, Bell and Jackiw [3] 



/IV 



£ pv pa^Xp^la • 



(147) 
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With this expression for f^, S(p 2 ) still has the form ( |142j ) but it will be quadrat- 
ically divergent and 



a 



71 



f=fa+ ~ A'-oo (148) 



where A is the ultraviolet cutoff. This means that Goldstone particles exist in the 
anomalous case but they are decoupled from any physical state. At the same time 
the Ward identity still makes sense because the product gf does not depend on 
/. Nevertheless, the concrete form of the Ward identity will change. The reason 
for this is again the Adler-Bell-Jackiw anomaly [3]. For the triangle diagram 



K> = / \ (149) 





the replacement of the pseudovector coupling by the pseudoscalar coupling gives 
an incorrect result: instead of the correct expression 



2rwy 5 a 

Pn&L = / \ + —epasykafap (150) 

which was obtained in [3], we get just the first term. We can try to write the 
Ward identity using (|150| ), but this turns out not to be necessary. The reason 
is that in this approach the axial current T^ a between gluonic states, defined 
symbolically by the relation 



^ = r? CT + > (i5i) 

h 9pa 



(where the term T^ a is defined diagrammatically and g pcr = y-r^ CT ) is just the 
transverse part of T^ a . 

rp 1 . = K - Vj f-K- ( 152 ) 

For the axial current between quark states T M can be written as 

T^ = F(q 2 ,q 1 ) + T^q 2 ,q 1 ). (153) 
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Here 



(154) 



is the same set of diagrams as in the non-singlet case and is the " axial current 
of gluons" 



9 



rp 

6 
/ \ 



+ 



rp 

6 

1 1 \ 



+ ■ 



(155) 



In the same way the Goldstone boson - quark interaction can also be divided 
into two parts. In the first part we can replace the pseudovector coupling by the 
pseudoscalar coupling. The second part is the longitudinal part of T^. Conse- 
quently, 




(156) 



and the Ward identity can be written in the form 



p,T - l fg = 1 ,G-\q 1 ) + G-\q 2 ) l5 . 



(157) 



This expression enables us to answer the question, what happens with particles 
like rj which were Goldstone states if we would not take into account that they 
can decay into two gluons. Let us consider the contribution of the massive pseu- 
doscalar flavour singlet particle r/ to the Ward identity ( |157| ). PpT^ has a pole 
corresponding to rj': 




9 V > i 

= V frf —n o ^ • 

j[T — P Z 



The Yukawa coupling to the Goldstone state has also a pole: 

{75, G- 1 } 



-fg 



n — p 



9rf 



158) 



(159) 
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The right-hand side of ( |157| ) has, however, no poles. This condition can be 
satisfied if 

» 2 fr,> = {j5,G- 1 (q)}<ZZ>g ri '. (160) 
The same Ward identity (|157|) gives at p 2 = 

frfQrf = {ils,G-\q)} (161) 

from what it follows that \j? for rj' is equal 

2 {^G-\g)}<CZ>{^G-\ q )} 

Vrf = 72 • ( 162 ) 

Jrf 

This means that rf acquired a mass due to the transition into a Goldstone boson 
which itself is decoupled. It is interesting to notice that at relatively small /1 2 
when the comparison between Ward identities for different values of p 2 (p 2 = fi 2 , 
p 2 = 0) makes sense, (|162|) gives us the same result as the equation ( |121| ) (without 
the additional point-like term which is present in the 7r-meson case ( |133| )). Indeed, 
in ( 121]) we can write 






0. 



The first two terms are equal to p (as in ( |133| )) and therefore 




P = \ J ~ -g<ZZ> g\ P =o, (163) 



in agreement with ( |162[ ). We see that the subtraction term in the ir- meson self 
energy, reflecting the quasi-local structure of the pion, disappears in the case of 
rj'. In this sense rf is a normal bound state without a point-like structure. 

The approach we presented here for the resolution of the [/(l)-problem is tech- 
nically very close to the approach developed by Veneziano [4] but the underlying 
physics differs essentially. In Veneziano's approach big long-range fluctuations 
are responsible for the rf mass. In our approach rf is a normal qq bound state 
with no local structure which would be responsible for its small mass if it were a 
Goldstone state. This local structure is destroyed by the decay on hard gluons. 



8 QCD with massive quarks 

We have discussed in detail an asymptotically free theory with massless fermions. 
We came to the conclusion that in order to obtain the correct spectrum of 
Goldstone particles, axial current conservation has to be imposed on the the- 
ory. QCD, however, contains massive quarks and the same spectrum of massive 
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quasi- Goldstone particles (the pseudo-scalar octet) as the theory with massless 
quarks. The problem is, how to impose the condition of axial current conservation 
on a theory which obviously does not conserve the axial current. 

In general, I don't know how to do this. For our real world, however, there is 
a natural possibility to solve the problem. 

In the real world QCD is part of the standard model describing strong, electro- 
magnetic and weak interactions. In the standard model all particles are supposed 
to be intrinsically massless and their masses appear as the result of symmetry 
breaking due to some kind of Higgs mechanism with or without elementary Higgs 
particles. The possibility of such a mechanism is guaranteed by the conservation 
of the left-handed SU(2) current j". For the matrix element of this current 
between any two quarks with momenta q 2 , qi we have the Ward identity 

P,K(q 2 , q x ) = \r a \{l - l5 )G-\ qi ) - G-\q 2 ) l -r a l -{l + 7s ). (164) 

In the case of massive fermions p^Ta contains the contribution of three Goldstone 
bosons responsible for the masses of W and Z° 

p,K = p,K - fg. (165) 

For large q 2 values G _1 contains a massive term Z~ 1 itiq. Hence, at large q 2 we 
have | | 

(fg)o = -^{t 75, moZ- 1 } + -[r a ,m Z _1 ]. (166) 

At small q 2 around the QCD scale fg is defined by the total quark mass m(q) 
which is for the light quark much larger than ttlq. Using the relations ( |164| ) - ( |166j ) 



we can calculate the the masses of the pseudoscalar octet in the same way as we 
did for T] f . Suppose there is a bound state of light quarks which is a pseudoscalar 
particle (the 7r-meson) with a finite mass fi. The pole corresponding to this 
particle will contribute to both terms in ( |165| ). With this pole contribution, p^T 1 ^ 
and fg can be written in the form 

P^a = < + (P^a)n,P = + (p^ a )n, P , (167) 

fl Z — p Z 



fy= }' h +(/j) v = ®^ + (/jW ( 168 ) 

pi — p 
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As in the case of rj', the condition for the pole cancellation gives us the value of 

n 2 = , 

At p 2 = we have 



fi 2 = Ufg)o<^>9,. (169) 



Ug* + (fg)n, P = ^ { r ^s, rh{q)Z x } + - [r a , rhZ x ] . (170) 

The Yukawa coupling g n of the bound state has to decrease at large q 2 . Therefore 
we have to identify (fg) n , P with (fg)o- Then 

Ugn = ]^T a -i^m s (q)Z- 1 (171) 

where m s (q) is the isotopically invariant part of the quark mass term behaving 
at q — > 00 as v 3 jq 2 . Hence, (fg)o ~ mo is defined by weak interactions with 
Goldstone states, and f n g n ~ m s is determined by pion exchange. The conclusion 
of these considerations is that in the case of massive quarks the conservation of 
the left-handed SU(2) current can play the same role for the calculation of the 
coupling and the masses of the flavour non-singlet pseudoscalar particles as in 
the massless case the conservation of the axial current. The result is that the 
masses become different from zero while the Yukawa coupling remains the same 
at least for small mo values. 

The expression ( |169| ) for fi 2 can be also obtained without discussing the cur- 
rent conservation, just by calculating S(p) as it is defined by ( |133| ) for the quark 



mass m s + m , keeping only the term linear in rriQ. It is in agreement with the old 
idea that mo influences the mass of the Goldstone state but not its wave function. 

Unfortunately, if we calculate /i 2 by using fll69| ) or (|133|) , we get a logarithmic 
divergence. For the 7r-meson mass the main divergent part equals 

ml = T ^Miyid) (172) 



4vr 2 f 2 Jul q 2 

(here v Q is of the order of Xqcd)- 

In the standard model the behaviour of m and z/ 3 can be calculated in a 
fantastically large region of q 2 : from 1 GeV up to the scale where one of the 
couplings of the standard model (g± - the U (l)-coupling, h - the Yukawa-coupling 
or A - the coupling of the self-interaction of the Higgs particles) becomes of the 
order of unity. In the reasonable case when the [/(l)-coupling gi is the first to 
become unity this scale is equal to A = 10 38 GeV. 

At q 2 > A 2 the behaviour of mo and u 3 is unknown. Because of this, m n is 
not calculable in principle and it has to be considered as an arbitrary parameter. 
If we, however, assume that at q 2 larger than A mo{q) and f 3 {q) become equal 
to zero, we will be able to calculate m n and, surprisingly, this calculation gives a 
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reasonable value for m n with A m 10 38 GeV [2]. The expression Q172 ) for m 2 is 
also in agreement with the naive expression 

2 2m - 

J IT 

with the important difference that now (\E , \[>) is defined not only by strong but 
also by weak interactions and it goes to infinity if the weak interaction is removed. 



9 The pion contribution to the equation for 
light quark Green's functions 

From the previous discussion it is obvious that the small mass pion contribution 
has to be included in the equation for the light quark Green's function. Fortu- 
nately this is very easy to do. Having in mind that now the diagrams contain not 
only the gluon contribution but also the emission and the absorption of pions we 
will find that as before, the main contribution to d 2 G~ 1 comes from the simplest 
diagram r rP ^ L , with the coupling {275, G^ 1 } at zero momentum instead of the 
gluon coupling d„G~ l . It leads to the following equation for the Green's function: 



d 2 G-\q) = gd,G- 1 (q)G(q)d,G- l (q) - {ij 5 , G~ 1 }G(q){ij 5 , G^ 1 } 



(173) 



The equation for bound states (|27|) has also to be changed. The correction comes 
from the diagrams 






Instead of (p7| ) we will have 
d 2 (j)(p } q) = 



g(q){A u (q 2 )d u (j)(p, q) + d v <j>{p, q)A v { qi ) - A v {q 2 )<j>{p, q)A v {q 1 )} + 



+ 



4vr 2 / 2 

J 7T 



{ l ^G-\q 2 )}G(q 2 )^cl>^G(q 1 ){ ll5 ,G- 1 (q 1 )} 



XG(q 1 ){i l5 ,G- 1 (q 1 )}^ - { ll5 ,G-\q 2 )}G(q 2 )^\ 



(174) 



Here A is the emission amplitude of the zero momentum pion in the transition 
of the bound state to the gg-pair. This amplitude has to be defined by the 
axial current conservation. There is another important quantity in this equation, 
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namely: /„■. In section |6] we have written an explicit expression for f 2 including 
only the gauge field contribution and ignoring the pion contribution. Now we 
include the pion contribution in the equation for the Green's function and, to be 
self-consistent, we have to do the same for f 2 . I was not able to carry out this in 
any order of the Yukawa coupling, but in the first order of g 2 the equation ( |137| ) 
is correct if one adds the diagrams 



^G- 1 ^G- 1 ditG -i dfiG -i 



{75, cr 1 }<^w^> {75, g- 1 } + {75, G ~ 1 ><C3^ > {75 ' G ~ 1} 



(175) 



As we see, the gluonic correction of the order of - and the pionic correction of the 
order of g 2 have the same diagrammatic structure. In order to estimate the value 
of these corrections, let us take just the contributions of zero momentum gluons 
and pions to them. It can be shown that the contribution of zero momentum glu- 



ons cancels in the last two diagrams of Ql37| ). Zero momentum pion contribution 
comes only from the first diagram of ( |1 75|) . Transferring the differentiation from 
the fermionic line to the pionic line in this diagram we obtain the contribution of 
zero momentum pions in the form 

{75, G- 1 }<^jT {75, G- 1 }-^ . (176) 
- Stt z 

The expression for f 2 which includes the zero momentum pion contribution is 
8/ - = /^p-rr{7 5 ,G- 1 }G{75,G- 1 }G3,A3 ( ,A + 

+ 84/(^ Tr({75 ' G " }G)1 - <177) 

It gives us an understanding of the scale of possible pion contributions. It is in- 
teresting to note, that (|177|) is not an expression in terms of the Green's functions 
but an equation for /„.. 

In the next paper we will show that the pion contribution changes essentially 
the structure of the equation for the Green's function. The new equation has a 
solution corresponding to the confined quark. At the same time the symmetry 
breaking solution will not necessarily survive (at least if g(0) is large). 
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